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Abstract 


I A theoretical analysis is conducted to study the effect of a tangential wall jet on the 

!. control of two dimensional separated flow past a circular cylinder. For a tangential 

■ wall jet, the mathematical model derived previously 1,2 is used and the vortex cloud 

[ method is adopted for the calculation of the external flow field. For certain limiting 

I 

j cases, the governing equations are simplified and closed forms of solutions for the 

! wall jet parameters can be obtained. 

It is observed that the wall jet is very efficient in reducing drag by delaying the 
L separation point in the range of small blowing strength. The suction force induced 

f 

I by the wall jet is negligible compared to the drag due to the external stream. 

However this suction force increases the drag when the blowing strength is large. 
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Nomenclature 


b half width of the wall jet 

c radius of cylinder 

C f friction coefficient 

Cfx jet blowing coefficient 

d horizontal distance of separation point 

G shape function 

K constant determined experimentally 

k constant determined from the velocity profile 

m number of bound vorticity panels 

N = kn 

n constant determined from the velocity profile 

p pressure 

q nondimensional variable for jet flow field 
R radius of curvature of body 

s, y surface coordinate system 

t nondimensional variable for jet flow field, time coordinate 
U free stream velocity 

u tangential component of velocity 

v normal component of velocity 

T vortex strength 

iv 


i- 


length scale of excess mass 
e eddy viscosity 

6 length scale of excess momentum 

p density 

<r complex coordinate 

C vortex core radius 

r viscous stress 

<f> velocity potential 

$ stream function 

0, vorticity distribution 

u> bound vorticity intensity 

< > time averaged quantity 

Subscripts 

0 zero pressure gradient 

e external stream 

j jet 

w wall 

m jet maximum velocity 


1. Introduction 


The wall jet may be used for the control of circulation around an airfoil. Here its use 
is investigated for drag reduction by delaying the flow separation in the presence of 
an adverse pressure gradient. Previous studies of the wall jet have not provided 
satisfactory results regarding the relationship between the blowing strength and the 
amount of drag reduction due to delayed separation. 

Recently a simple method * has been derived for the two dimensional turbulent 
wall jet for an incompressible flow, including the effects of the external stream. In 
this study the wall jet model is applied to the case of the flow around a circular 
cylinder, which can be easily transformed to other body shapes. The external 
stream which is imposed on the jet model is calculated using the vortex cloud 
method. 6,7,8 As a preliminary study for the wall jet application, only the mean 
quantities averaged for a long time are calculated. Typically attention is focused on 
the effects of the wall jet strength on the variation of drag force and the displacement 
of the separation point. 


2. Mathematical Model 


The mathematical model for a 2 dimensional incompressible turbulent wall jet was 
derived in Ref. 1 and 2. In this section the model for the wall jet is briefly reviewed 
and the governing equations for the wall jet are modified for the limiting case of jet 
boundary conditions. The vortex cloud method is employed to calculate the exter- 
nal flow field. Due to the inviscid nature of the vortex cloud method appropriate 
assumptions are introduced for the calculation of surface pressure distribution. 

2.1 Tangential Wall Jet 

2.1.1 Mean Velocity Profiles 

The mean velocity profile of the wall jet is approximated as comprising two parts; 
(1) an inner flow adjacent to the curved surface having the characteristics of a 
turbulent boundary layer, (2) an outer flow which merges with the external stream 
and is typical of a free turbulent jet. These profiles are shown in Fig. 1. These two 
parts can be written as follows 

“ = um [ 2 (£) - (£) _ for y - ym (211) 

u ~ = (um - «e) sech 2 k [7 — — for y > y m (2.1.2) 

L 0 — l lm . 

where n depends on the Reynolds number of the jet, and lies in the range 5 to 9. 
The definition of 6, half width of the jet, determines the constant k as 

k = tanh -1 4= = 0.8814 (2.1.3) 

V2 


2 



The location of the maximum velocity y m can be found by matching the second 
derivative of the velocity profiles in the two regions. This results in 


Vm 

b 



(2.1.4) 


In this study it is assumed that the external velocity u e is smaller than the jet 
maximum velocity u m throughout the jet flow field. 

2.1.2 Governing Equations for the Wall Jet 


With the assumptions that 

(a) the contribution to the overall momentum balance of the shear stress at the 
wall is small 


(b) the effects of the wall curvature are small except where it influences the gra- 
dients of pressure and shear stress, 

the equations of continuity and momentum for a 2-D incompressible flow are 


du dv 
ds + dy 


= 0 


(2.1.5) 


du dv d { p\ d ( t\ 

U ds dy ds \p/ dy VpJ 

where the pressure and the shear stress are determined from 

— ft) m2 

dy\p) R + y 


( 2 . 1 . 6 ) 


(2.1.7a), 


r 

P 



(2.1.76) 


where R is the radius of curvature of the body surface 


and e is an eddy viscosity. 
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The equation obtained by combining Eqs. (2.1.5) and (2.1.6) can be written in 
integral form as 


d ( Q 2\ i d u e /<-* \ T w 

5ll Su d + 17 (ii e) = “T 


( 2 . 1 . 8 ) 


where 9 and 6* are the excess of the momentum and mass thicknesses in the wall 
jet respectively, and are defined as follows: 


-fife-)* 


(2.1.9a) 


(2.1.96) 


Substitution of the velocity profiles given in Eqs. (2.1.1) and (2.1.2), ignoring higher 
order terms and the small shear stress at the wall results in 


d_ 

ds 


b ( Um — We) 


^w m + — 


+ 


— 6(fc + 1) (um — We) 


du e 

ds 


= 0 ( 2 . 1 . 10 ) 


Another equation can be obtained from the momentum equation at y = y m , with 
an appropriate eddy viscosity model, as 


1 d { 
Kds V U 


2 

m 



( w m -Ue f 
+ b 


3 1 

/. ym\~ l n(w m - w e )2 w^ _ 
\ b ) + k R 


( 2 . 1 . 11 ) 


where K is a constant to be determined experimentally. 

These two equations, (2.1.10) and (2.1.11) are the governing equation for the two 
jet parameters 6 and u m . 


2.1.3 Curvature-Induced Suction at the Wall 


Even though the wall jet reduces drag by delaying separation, it causes a rearward 
suction force due to attachment of the jet to the surface before separation - the 
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Coanda effect. This amount of suction force can be calculated by integrating the 
radial momentum equation given by Eq. (2.1.7a) as 


Ps 


r°°u 

~ Pe = ~ / ~ 

JO 


°°’' 2 - u e 2 


dy 


i o R + y 

where s represents the value at the wall surface. 

With the definition of 9 and 8* it can be reduced to the approximate result 


( 2 . 1 . 12 ) 


Ps ~ Pe — 


R+ 2 


(2.1.13) 


2.1.4 Wall Jet Separation 


For a turbulent boundary layer the separation condition 1 may be written as 


»Hi) 


T w,0 


= 5 


(2.1.14) 


Using the definition of 9 and the velocity profiles given in Eqs. (2.1.1) and (2.1.2), 
this condition for the wall jet with large n can be written as follows 


LA. ( p \ = 

Kds \p) e 2 b K 
where Cf is defined as 


1 + kn 


\ u m) 


i- 

-1 


l_d_ 

+ Kds 


(» + S') ug 
R + 2 


C f = T 


Tw 


yP u m 


(2.1.15) 


(2.1.16) 


2.1.5 Dimensionless Forms of the Equations 


The following nondimensional variables are used 

8 (u m — u e )u m 


q = 


c u 


oo 


, “e _ 

t = , It 

Um 


Ue 

Uoo 


(2.1.17) 
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where c is cylinder radius and u is known a priori from the external flow field. 


In the case of a cylinder the radius of curvature R = c. With these variables 
and normalizing all the length scales with the cylinder radius c, Eqs. (2.1.10) and 
(2.1.11) become 



From now on and in the equations shown above, even though all the symbols rep- 
resenting the length scale are quantities normalized with c, the same symbols are 
used for convenience. 


These two equations, (2.1.18) and (2.1.19), are solved subjected to the conditions 

6(0) ( Uj - Ue(0)) Uj 


?U=0 = 


U 


oo 


£|s=0 — 


“e(0) 

Uj 


(2.1.20 a) 
(2.1.206) 


where uj is the velocity at the jet slot. 

When 6(0) becomes 0, the ideal case of a jet slot, uj becomes infinite to maintain 

a finite momentum injection and these boundary conditions reduce to 

6(0)u| _ 


,(0) = — i-L = C„ 


u 


oo 


<(o) = ~ o 




(2.1.21a) 


(2.1.216) 


Rearranging Eqs. (2.1.18) and (2.1.19), the following set of ordinary differential 
equations can be obtained 

* ‘ (l _ *) + + ()f ( 2.1.22) 

ds 2 it 2 V / 2 qt N (1 — t) 2 2A: 2 
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(2.1.23) 


dq _ qt (1 - t 2 ) (fir 2 _ ir (l-<) 3 fr 2 l + N (1-t)? 
ds 5u 2 1 + 0.4< ds 5(1 + 0.4t) ^ 

gtifiV (1 - <)f 3(Jb + l) 1 du 2 

5fc 2 1 + 0.4t ” 8 1 + 0.4t ds 


where N = kn 

These two equations are solved simultaneously using the 4^ order Runge-Kutta 
method. 

Eq. (2.1.13) for the pressure at the wall can be written in terms of q and t 


Ps ~ Pe 

\pUoo 2 


3k q 


1 + (I + f + !)) i 

1 1 1 t2q 

1 + 2i*(ir^ 


(2.1.24) 


and the separation condition given by Eq. (2.1.15) is 


1 1 (dp\ 5 /*C/\ ^(l -t) r ( ai 

A-lp Uoo 2 UJ ( 2k \ K ) t*q l 1+ ( H 


L—— 

+ K3kdl 


9 (* + (I + I (* + 1 )) *) 


1 1 1 t2q 


(2.1.25) 


2.1.6 Reduced Equation around the Singular Point 


The above two equations (2.1.22) and (2.1.23) show a singular behavior around the 
jet exit where the boundary conditions given by Eqs. (2. 1.21a, b) are imposed. To 
avoid this singularity, the two equations can be expressed around the region of small 
t as 


dq ~ K l + N u 2 
ds ~ 5 N t 


(2.1.26) 
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(2.1.27) 


dt_ ~ K l + N u 2 
ds ~ 2 N qt 

Combining theses two equations with the condition given by Eq. (2.1.21) the closed 
form solutions for q and t can be obtained as 

q — Cfx exp(— 0.4t) (2.1.28) 

and 

(2X29) 

These* solutions are applied to the region around the jet exit. 

Eq. (2.1.29) shows that t increases rapidly, which means that u m decreases 
quickly. This indicates that in the vicinity of the slot, the escaped jet diffuses 
rapidly. This phenomena causes a sudden increase of b, as can be seen in Eq. 
(2.1.28). 

2.2 Calculation of the External Flow Field 

In this study, it is assumed that the wall jet affects the external flow field by changing 
the separation position and, thereby, the external pressure distribution. Thus for 
an assumed separation point the external flow field can be calculated separately 
from the jet flow field. The external flow field is calculated using the vortex cloud 
method, such that at each time step a new vortex is generated from the separation 
point. 

2.2.1 Vortex Cloud Method 

It is well known that the vortex cloud method is very efficient for the calculation of 
a two dimensional unsteady separated flow field. 
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For a 2-D flow field, the governing equation and the boundary conditions, in 
terms of the stream function can be written as 0 

V 2 $ = -ft ' (2.2.1) 

^ = const. on surface (2.2.2a) 

\t r = U S<r in the far field (2.2.2 b) 

where 5 represents the imaginary part of a complex value and ft is the vorticity 
distribution in the flow field which is replaced by a cloud of vortices of finite core 
size to avoid infinite induced velocities. 

Introducing the core radius <; and the vortex strength T,- for the i th vortex, ft can 
be written as 

< 2 - 2 - 3 > 

where G is a shape function common to all vortices. 

In this study the core radius <; is kept constant throughout the calculation proce- 
dure. 

The shape function and the associated velocity distribution used in this study 
were derived by Spalart® 

^ _J° 0>1 

* G ~ { z(e 2 - 1) 2 e < i 

2ttuc _ f \ Q > 1 

r, “ \ e (z - 3<? 2 + e 4 ) e<i 

where g = |<7 — <r,-| /<;. 

A new discrete vortex is generated from the prescribed separation point at each 
time step such that the flow field including the new vortex satisfies an appropriate 


(2.2.4) 

(2.2.5) 
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Kutta condition at the separation point. In an unsteady flow field the vorticity is 
continuously supplied to the surrounding fluid from a separation point such that 
the strength and the position of the shed vorticity allow for pressure continuity to 
exist. In this study a simple Kutta condition 6 is used to reduce computational cost. 
Using this technique the new vortex position and strength are determined as 

r n ew = 2u s 2 At (2.2.6) 

Anew = & 8 d" CgUsAt (2.2.7) 

where s represent the quantity at the separation point and e s is the tangential 
unit vector at the separation point. At is the time increment of the calculation 
procedure 

2.2.2 Numerical Implementation 


The cylinder surface is represented by vorticity panels of linearly varying strength. 
If the surface is divided into m vorticity panels and there are n vortices in the 
separated flow field, as shown in Fig. 2, the stream function at the control point i 
is 


*i = j> u(°j)K(<Ti,(Tj)daj 
1 " 

+ lo S I* 7 *' - a k\ + USv, i 


( 2 . 2 . 8 ) 


Jk=l 


where u{<jj) is the bound vorticity density at the panel and K(<r^<rj) are the 
geometric influence coefficients, which depend only on body geometry. 

The first term on the right-hand side of Eq. (2.2.8) is the contribution of the 
bound vorticity panels, and the second is the contribution of the discrete vortices 
representing the separated flow field. The last term is due to the free-stream. 
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Boundary condition (2.2.2a) gives 


= ^2 = • • • = = const. = c 

Conservation of vorticity requires 



(2.2.9) 


( 2 . 2 . 10 ) 


There are m + 1 unknowns; m u > ’s and c, the value of the stream function at the 
wing trace. Eqs. (2.2.8), (2.2.9) and (2.2.10) constitute a linear system of m + 1 
equations. 

The calculation procedure by solving the equations (2.2.8), (2.2.9) and (2.2.10) 
is repeated for a time period large enough to obtain mean flow field quantities. 
At each time step a new vortex is generated from the prescribed separation point 
and convected with the velocity induced by the neighboring vortices and the bound 
vorticity panel. 

For computational efficiency the following techniques are used to reduce the num- 
ber of vortices in the flow field and to avoid numerical singularities. 

(1) if the vortices come closer to each other than a given value the two vortices are 
allowed to merge into one and the position and strength Eire determined such 
that linear and angular momentum before and after merging are conserved. 

(2) if the vortices approach closer to the surface than a given value the vortices 
are made to disappear 

(3) if the velocity induced on the surface by the vortices located far downstream 
is less than a given value the vortices are made to disappear 


2.2.3 Calculation of Mean Pressure and Velocity 
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From the unsteady Bernoulli equation, the pressure coefficient can be written as 



( 2 . 2 . 11 ) 


The velocity at the surface can be obtained from the definition of bound vorticity 
intensity as 7 


u e _ u> 

a ~u 


( 2 . 2 . 12 ) 


The mean quantities axe obtained as follows 


<C D > 




C p dt 


l fT u 

< “ >= f jL vz dt 


where T is the toteil time during the calculation proceeded. 


(2.2.12a) 

(2.2.126) 


2.3 Overall Calculation Procedure 

Once the mean quantities of the external streeim axe obtained from the vortex 
cloud method, these values are substituted into equations (2.1.22) and (2.1.23) to 
calculate the jet parameters g, t and Cp , which satisfy the separation condition 
(2.1.25) at the prescribed separation point. Since the boundary condition (2.1.21) 
causes a singularity at the jet exit, the approximated solutions (2.1.28) and (2.1.29) 
for q and t are applied to the first bound vorticity panel from the jet exit. The 
parameters, q and t, are calculated from the jet exit to the separation point. 

Drag can be calculated by integration of the surface pressure, which can be 
obtained from Eqs. (2.1.24) and (2.2.11). The vortex cloud method provides a 
realistic pressure distribution when the vortices are located far from the surface. 6 
For the present problem, however, the separated vortices stay close to the surface 
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in the aft-separation region, where the viscous effect are also significant.® Here the 
pressure is assumed constant after separation, an assumption that is adopted in 
many studies for the cylinder problem. 
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3 Results and Discussion 


In this study it is assumed that the wall jet affects the external flow field by changing 
the separation point. Although this assumption is valid only for a small blowing 
strength less than 0.05, the blowing strength is allowed to vary to larger values 
to investigate the quantitative trends of the flow quantities. The value of n, the 
parameter for the mean velocity profile, is taken as 7 and Cf , the friction coefficient, 
is chosen as 0.01. The overall flow pattern and jet slot position are shown in Fig. 
3. 

3.1 Results 

Fig. 4 shows the wake patterns behind a cylinder for four different separation 
positions. From the figure, it can be seen that the wake width reduces, while the 
frequency increases as the separation point is delayed. In the case of d/c = 0.26, 
(which is the natural separation observed experimentally without jet blowing) the 
separation angle is 106 deg. and only three vorticity clusters exist in the given 
downstream region, but for d/c = 0.8 four of those exist. The variation in width of 
the wake is also clearly visible from the figures. 

Fig. 5 shows the amount of suction force due to the wall jet, given by Eq. 
(2.1.24). The pressure difference across the jet is very small. The blowing strength 
for this case is very large, its value is almost 1, since for small values of blowing 
strength the pressure difference can be barely observed. 

Fig. 6 shows the variation of pressure distribution on the cylinder surface for 
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different separation points. As the separation point is delayed, the pressure distri- 
bution is recovered to that of the potential calculation without separation except 
for the small amount induced by the wall jet. 

Fig. 7 shows the velocity distribution in the wake at a downstream location. One 
is for the case of natural separation, d/c = 0.26, and the other is for djc — 0.7. 
This figure clearly shows the difference in the momentum deficit of the wake, which 
causes a drag on the cylinder. The momentum deficit is reduced significantly with 
a delay of the separation position. 

Fig. 8 shows the variation of jet parameters for different blowing strengths. The 
limit solutions given as Eqs. (2.1.28) and (2.1.29) are applied only to the first panels 
after the jet exit and the full equations, (2.1.22) and (2.1.23), are applied to the 
rest of the region. Fig. 8a shows the variation of t , the ratio of external velocity to 
jet maximum velocity, t increases rapidly close to the jet exit, indicating a sudden 
decrease of the jet maximum velocity due to diffusion of the jet. The diffusion in 
the jet also can be seen in Fig. 8b. Fig. 8a also shows that far from the jet slot 
the maximum velocity relative to the external velocity reduces with a constant rate 
regardless of blowing strength, while q varies with a different slope for different 
blowing strengths as can be seen in Fig. 8b. 

Fig. 9 shows the variation of drag coefficient with blowing strength. In the region 
of small blowing strength, a slight change of Cp causes a significant decrement of 
drag. For instance when = 0.06, which corresponds to d/c ~ 0.6, the drag is 
reduced to almost 1/3 of that without blowing. This phenomena shows the very 
efficient mechanism of wall jet blowing with a small amount of momentum injection. 
For large blowing strength, the drag does not change much with further increments 
of Cfi . As can be seen in the figure, drag even increases with , due to the suction 
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force induced by the jet. Also in this region the assumption that the jet does not 
affect the external flow field could be violated. 

Fig. 10 shows the amount of separation delay with respect to Cp . This figure 
shows trends similar to the variation of Cp . Small increments of Cp provide large 
separation delays in the small Cp region. For the large Cp region, even with a large 
variation of Cp, the separation point does not change appreciably. 

3.2 Concluding Remarks 

An analytical study of a tangential wall jet applied to the control of separated flow 
past a cylinder has been conducted. The external flow field has been calculated 
using the vortex cloud method and a mathematical model for a 2-D incompressible 
turbulent wall jet was applied inside the wall jet flow field. 

The vortex cloud method provides a good prediction of the external flow quanti- 
ties, such as surface velocity distribution and wake pattern. As the separation point 
is delayed, the width of the wake is reduced and the frequency of the Karman- vortex- 
street increases. 

It is observed that the thin tangential wall jet is very efficient in drag reduction 
by delaying separation for small blowing strengths, but is not efficient for large 
blowing strength. For large blowing strength drag even increases due to suction of 
the jet. 

Although this study provides a reasonable indication of the relation between Cp 
and Cq, it is desirable to have experimental data to validate the analysis. In this 
study only mean quantities of the external stream are imposed on the jet model to 
separate the calculation process between the external stream and the jet flow field. 
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The wall jet model is limited by the assumption u m > u e throughout the jet 
flow field. Whereas experiments 4,11 have shown that u e may exceed u m due to 
severe adverse pressure gradients around the blunt trailing edge. The assumption, 
u m > , can be eliminated, however, by the introduction of a additional parameter 

in the mean velocity profile as in Ref. 3. 

This study also assumes that the jet does not affect the external flow field ex- 
cept by changing the separation point. This assumption is valid for small blowing 
strength; for large blowing strength, the external flow could be influenced by the 
jet as suggested in Fig. 11. 11 
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Fig. 1 Velocity profile around jet flow field 
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Fig. 2 Bound vorticity panels and separated vortices 
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Fig. 7 Momentum deficit at x/c » 2 for two different separation points 
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Fig. 8 a Variation of t from the jet slot to the separation point 



Fig. 8 b Variation of q from the jet slot to the separation point 
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